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Abstract 

In this article, we study the ^ + doubly heavy baryon states £Iqq and Sqq by sub- 
tracting the contributions from the corresponding i doubly heavy baryon states with 
' the QCD sum rules, and make reasonable predictions for their masses. Those doubly 

heavy baryon states maybe observed at the Tevatron, the LHCb and the PANDA. 
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■ 1 Introduction 

In 2002, the SELEX collaboration reported the first observation of a signal for the doubly 
charmed baryon state in the charged decay mode — > A+K~ir + [lj, and confirmed 
later by the same collaboration in the decay mode — > pD + K~ with the measured 



mass Ms = (3518.9 ± 0.9) MeV [2]. However, the Babar and Belle collaborations have 
J> ■ not observed any evidence for the doubly charmed baryons in e + e~ annihilations OH]. 

The charmed and bottom baryons which contain one (two) heavy quark(s) are particularly 
interesting for studying dynamics of the light quarks in the presence of the heavy quark(s), 
Tj" ■ and serve as an excellent ground for testing predictions of the quark models and heavy 

quark symmetry. There have been several approaches to deal with the doubly heavy 
baryon masses, such as the relativistic quark model 0[6], the non-relativistic quark model 
[3 IH El [10] , the three-body Faddeev method [UJ, the potential approach combined with 
the QCD sum rules |12| . the quark model with AdS/QCD inspired potential [13], the MIT 
bag model [13], the full QCD sum rules |15[ll6j. and the Feynman-Hellmann theorem and 
^ ' semiempirical mass formulas [T7], etc. 

The QCD sum rules is a powerful theoretical tool in studying the ground state heavy 
baryons [18l QjJJ [20]. In the QCD sum rules, the operator product expansion is used 
to expand the time-ordered currents into a series of quark and gluon condensates which 
parameterize the long distance properties of the QCD vacuum. Based on the quark- 
hadron duality, we can obtain copious information about the hadronic parameters at the 
phenomenological side [181 [19] . There have been several works on the masses of the heavy 
baryon states with the full QCD sum rules and the QCD sum rules in the heavy quark 
effective theory (one can consult Ref.|21] for more literatures). 

In Ref. [22], Jido et al introduce a novel approach based on the QCD sum rules to sep- 
arate the contributions of the negative parity light flavor baryon states from the positive 
parity light flavor baryon states, as the interpolating currents may have non-vanishing 
couplings to both the negative parity and positive parity baryon states [23]. Before the 
work of Jido et al, Bagan et al take the infinite mass limit for the heavy quarks to sep- 
arate the contributions of the positive parity and negative parity heavy baryon states 
unambiguously [24] . 
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In Refs. [25 | [26 | [2T] . we study the heavy baryon states Qq, E'q, Sq, Qq, Eq and £g with 

the full QCD sum rules, and observe that the pole residues of the | heavy baryon states 
from the sum rules with different tensor structures are consistent with each other, while 
the pole residues of the ^ + heavy baryon states from the sum rules with different tensor 
structures differ from each other greatly. In Refs. [211 EH]> we follow Ref . |22j and study 
the masses and pole residues of the ^ heavy baryon states Qq, E'q, Eg, Ag and Eq by 
subtracting the contributions of the negative parity heavy baryon states to overcome the 
embarrassment. Those pole residues are important parameters in studying the radiative 
decays £l* Q -> Q Q -f, E* Q -)■ E' Q j and Eg -> £q7 [27J [29], etc. 

In this article, we extend our previous works to study the ^ + doubly heavy baryon 
states Eqq and Qqq with the QCD sum rules. 

The article is arranged as follows: we derive the QCD sum rules for the masses and 
the pole residues of the doubly heavy baryon states Hgg and Qqq in Sec. 2; in Sec. 3, we 
present the numerical results and discussions; and Sec. 4 is reserved for our conclusions. 



2 QCD sum rules for the Eqq and Q,qq 

The ^ + doubly heavy baryon states Eqq and £Iqq can be interpolated by the following 
currents Jz.(x) and Jq{x) respectively, 

J B {x) = e^Qf^C^Qjixh^qkix), 

J n (x) = e^Qf^C^Q^x^Skix), (1) 

where the Q represents the heavy quarks c and b, the i, j and k are color indexes, and the 
C is the charge conjunction matrix. In the heavy quark limit, the doubly heavy baryon 
states can be described by the (heavy) diquark- (light) quark model [12]. 

The corresponding | doubly heavy baryon states can be interpolated by the currents 
J_ = ^7 5 J + because multiplying 275 to J + changes the parity of J + [22], where the J + 
denotes the currents Jq(x) and Js(x). The correlation functions are defined by 

n±(p) = ijd 4 xe^ x (0\T{j±(x)J±(0)}\0), (2) 

and can be decomposed as 

U ± (p) =filx{p 2 ) ±U (p 2 ) , (3) 

due to Lorentz covariance. The currents J+ couple to both the positive parity and negative 
parity baryon states [23], i.e. (0| J + \B-)(B~\ J+|0) = -75 (0| J-\B~ ) (B~\ J_| 0)75, where 
the B~ denote the negative parity baryon states. 

We insert a complete set of intermediate baryon states with the same quantum numbers 
as the current operators J+(x) and J-{x) into the correlation functions n + (p) to obtain 
the hadronic representation [181 EH] ■ After isolating the pole terms of the lowest states of 
the doubly heavy baryons, we obtain the following result [22J: 
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where the M± are the masses of the lowest states with parity ± respectively, and the A± 
are the corresponding pole residues (or couplings). If we take p = 0, then 
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the -A(po) + B(po) and -A(po) — -B(Po) contain the contributions from the positive-parity 
and negative-parity baryon states respectively. 

We calculate the light quark parts of the correlation functions n + (p) in the coordinate 
space and use the momentum space expression for the heavy quark propagators, i.e. we 
take 
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where 



( s a ' ), then resort to the Fourier integral to transform the light 



quark parts into the momentum space in D dimensions, take p = 0, and use the dispersion 
relation to obtain the spectral densities p A (po) and p B (po) (which correspond to the tensor 
structures 70 and 1 respectively) at the level of quark-gluon degrees of freedom. Finally 



we introduce the weight functions exp 
rules, 
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the sq are the threshold parameters, T 2 are the Borel parameters, aj 
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With the simple replacements, (ss) — > (qq) and m s —¥ 0, we can obtain the corresponding 
spectral densities of the Hqq at the level of quark-gluon degrees of freedom. 



3 Numerical results and discussions 



The input parameters are taken as (qq) 



(qg s crGq) = m^(qq), (sg s aGs) 



rriQ(ss) 



-(0.24 ± 0.01 GeV) 3 , (ss) = (0.8 ± 0.2) (qq), 
-- (0.8±0.2)GeV 2 EHUD], (^§^) = (0.012± 
= (1.35±0.10) GeV and m b = (4.7±0.1) GeV 



has been updated from time to time, and 



0.004) GeV 4 [31], m s = (0.14±0.01) GeV, m c = 
[32] at the energy scale [i = 1 GeV. 

The value of the gluon condensate ( a s GG ) 
changes greatly [20J. At the present case, the gluon condensate makes tiny contribution 
(see Table 1), the updated value ( s ^) = (0.023 ± 0.003) GeV 4 [20] and the standard 
value (£*§£) = (0.012 ± 0.004) GeV 4 [31] lead to a slight difference and can be neglected 
safely. 

The Q-quark masses appearing in the perturbative terms are usually taken to be 
the pole masses in the QCD sum rules, while the choice of the m,Q in the leading-order 
coefficients of the higher-dimensional terms is arbitrary \20\ |33j. The MS mass m c (m 2 ) 



relates with the pole mass m c through the relation m c (m; 
In this article, we take the approximation m c (m 2 



m. 
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in, 



without the a, corrections for 



1.27 



+0.07 
-0.11 



GeV 



consistency. The value listed in the Particle Data Group is m c (m 2 ) 
it is reasonable to take rh c = m c (1 GeV 2 ) = (1.35 ± 0.10) GeV. For the b quark, the MS 
mass m(,(m 2 ) = 4.201q'q7 GeV [32] , the gap between the energy scale fx = 4.2 GeV and 
1 GeV is rather large, the approximation ~ ra b (rn^ ) m 6 (l GeV 2 ) seems rather crude. 
It would be better to understand the quark masses m c and m\, we take at the energy scale 
/x 2 = 1 GeV 2 as the effective quark masses (or just the mass parameters). 

In calculation, we neglect the contributions from the perturbative 0(a™) corrections. 
Those perturbative corrections can be taken into account in the leading logarithmic ap- 
proximations through the anomalous dimension factors. After the Borel transform, the 
effects of those corrections are to multiply each term on the operator product expansion 



side by the factor, a \^fi) ' where the Tj is the anomalous dimension of the in- 

terpolating current J(x), and the To n is the anomalous dimension of the local operator 
O n (0), which governs the evolution of the vacuum condensate (O n (0))« with the energy 
scale through the renormalization group equation. 

If the perturbative 0(a s ) corrections and the anomalous dimension factors are taken 
into account consistently, the spectral densities in the QCD side should be replaced with 
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where the A(p^ : m 2 ^) and -B(pg,mg) are some notations for the coefficients of the per- 
turbative corrections, the average virtuality of the quarks in the correlation functions is 
characterized by the Borel parameter T 2 . We cannot estimate the corrections and the un- 
certainties originate from the corrections with confidence without explicit calculations. In 
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Figure 1: The contributions of the pole terms with variations of the Borel parameters T 2 , 
the A, B, C and D correspond to the channels E cc , S\ c , H&& and f2&& respectively, the (5 
corresponds to the central values of the threshold parameters, the energy gap among a, (3 
and 7 is 0.1 GeV. 



the case of the correlation function for the proton, the perturbative 0(a s ) corrections can 
change the numerical values of the mass and the pole residue considerably and improve 
the agreement between the theoretical estimation and the experimental data [301 134j . In 
the present case, including the a s corrections maybe improve the predictions. 

In this article, we carry out the operator product expansion at the special energy 



scale /x 2 = 1 GeV 2 , and set the factor 



2r./-r 0n 

- Tjjn ~ 1 for consistency, as the a s 

corrections have not been calculated yet. Such an approximation maybe result in some 
scale dependence and weaken the prediction ability. In this article, we study the J p = ^ + 
doubly heavy baryon states systemically, the predictions are reasonable as we take the 
analogous criteria in those sum rules. 

In the conventional QCD sum rules [18U19] . there are two criteria (pole dominance and 
convergence of the operator product expansion) for choosing the Borel parameter T 2 and 
threshold parameter sq. We impose the two criteria on the doubly heavy baryon states to 
choose the Borel parameter T 2 and threshold parameter sq. 

In Fig.l, we plot the contributions from the pole terms with variations of the Borel 
parameters T 2 and the threshold parameters sq. The pole contributions are larger than 
(or about) 50% at the values which are denoted by the vertical lines for central values 
((3) of the threshold parameters sq. From the figure, we can set the upper bound of 
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the Borel parameters T^ ax , T^ ax = 3.8 GeV 2 , 4.0 GeV 2 , 9.1 GeV 2 and 9.3 GeV 2 in the 
channels H cc , f) cc , Hj& and f2^, respectively. The lower bound of the Borel parameters 
Tmin can be determined in the regions where the contributions from the perturbative 
terms are larger than (or equal) the ones from the quark condensates, T^ in = 2.8 GeV 2 , 
1.6 GeV 2 , 7.7 GeV 2 and < 6.5 GeV 2 in the channels E cc , Q cc , Hjj,& and f^, respectively. 
The convergent behaviors in the channels £Iqq are better than the corresponding ones in 
the channels Hqq, this is mainly due to the fact that the values of the quark condensates, 
\{qq)\ > \(ss)\- The Borel windows T^ in — in the channels Hqq and Qqq overlap 

with each other. In this article, we can take the uniform intervals for the Borel windows 
(IBW), i.e. IBW = 1.0 GeV 2 and 1.4 GeV 2 in the doubly charmed and doubly bottom 
channels respectively, which warrant the pole contributions are analogous. The values of 
the threshold parameters sq and the Borel parameters T 2 are shown in Table 1, from the 
table, we can see that the two criteria of the QCD sum rules are fully satisfied [18t [19]. 

In Ref.[2T], we study the \ + sextet heavy baryon states flf,, J7 C , E' b , E' c , £{, and S c 
by subtracting the contributions from the corresponding negative parity heavy baryon 
states with the QCD sum rules, the predicted masses are in good agreement with the 
experimental data for the well established mesons, f^, Q c , rl c and E c . In those sum 
rules, the contributions from the pole terms are about (45 — 65)% and (45 — 80)% for the 
bottom and charmed baryon states respectively. In this article, we take analogous pole 
contributions, see Table 1, the predictions are reasonable although the contributions from 
the perturbative continuum are somewhat large. 

In Fig. 2, we plot the predicted masses with variations of the threshold parameters sq. 
From the figure, we can see that the predicted masses are not sensitive to the threshold 
parameters, although they increase with the threshold parameters. In calculation, we take 
uniform uncertainties for the threshold parameters, 5 V /^ = ±0.1 GeV. 

From Table 1, we can see that the contributions from the quark condensates are large, 
even comparable with the perturbative terms, this is an indication of the non-perturbative 
origin of the masses of the baryon states. In the chiral limit, the spectral densities on the 
QCD side come from the quark condensates only. An astonishingly simple expression can 
be obtained in case of the proton [35j 136]. 

M p = ^/-87r 2 (^) M=1Gc v ~ lGeV. (12) 

Taking into account all uncertainties of the relevant parameters, we can obtain the val- 
ues of the masses and pole residues of the doubly heavy baryon states Hqq and Qqq, which 
are shown in Figs. 3-4 and Tables 2-3. In Table 2, we also present the predictions of other 
theoretical approaches and the values of the experimental data, the present predictions 
are consistent with them. 

The fractions 
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are shown explicitly in Fig.5. At the value T 2 = (2.5 - 4.0) GeV 2 , R = (-7 ~ 17)% in the 
doubly charmed baryon channels; and at the value T 2 = (7.5 — 9.5) GeV 2 , R = (—4 ~ 10)% 
in the doubly bottom baryon channels. The contributions from the negative parity doubly 
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T 2 (GeV 2 ) 




pole 


perturbative 


(??) 






2.8-3.8 


4.2 


(46 - 78)% 


(48 - 55)% 


(43 - 50)% 


»2% 




3.0-4.0 


4.3 


(45 - 75)% 


(67 - 72)% 


(26 - 30)% 


«2% 




7.7-9.1 


10.8 


(46 - 69)% 


(48 - 53)% 


(46 - 51)% 


< 1% 




7.9-9.3 


10.9 


(46 - 68)% 


(69 - 72)% 


(28 - 31)% 


< 1% 



Table 1: The Borel parameters T 2 and threshold parameters sq for the doubly heavy 
baryon states, the "pole" stands for the contribution from the pole term, and the "per- 
turbative" stands for the contribution from the perturbative term in the operator product 
expansion, etc. In calculating the contributions from the pole terms, we take into account 
the uniform uncertainties of the threshold parameters, S^j^ = ±0.1 GeV. 
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Table 2: The masses M(GeV) of the doubly heavy baryon states, where the star * 
denotes the masses of the negative parity doubly heavy baryon states predicted by the 
non-relativistic quark model. 
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0.311 ±0.077 



Table 3: The pole residues A+ of the doubly heavy baryon states. 



8 




Figure 2: The masses M of the doubly heavy baryon states with variations of the threshold 
parameters so, the A, B, C and D correspond to the channels S cc , Q cc , and 
respectively, the Borel parameters T 2 are taken to be the central values. 
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Figure 3: The masses M of the double heavy baryon states, the A, B, C and D correspond 
to the channels H cc , Q cc , 3^ and respectively. 
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Figure 4: The pole residues A + of the doubly heavy baryon states, the A, B, C and D 
correspond to the channels H cc , VL CC , and respectively. 
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Figure 5: The ratios between the contributions from the negative parity and positive 
parity doubly heavy baryon states with variations of the Borel parameters T 2 , the A, B, 
C and D correspond to the channels H cc , f2 cc , H&& and Q&6 respectively. 
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heavy baryon states are not very large, although they are considerable in some sense. In 
this article, the central values of the threshold parameters are sq = 4.2 GeV, 4.3 GeV, 
10.8 GeV and 10.9 GeV in the channels H cc , fi cc , and fi^ respectively, which are larger 
than the masses of the corresponding negative parity doubly heavy baryon states, see 
Table 2. If we take the tensor structures p and 1, the contributions from the negative 
parity baryon states are included in. In the case of the doubly bottom channels, the 
threshold parameters are slightly larger than the thresholds of the corresponding negative 
parity baryon states, the contaminations are small, see Fig. 5. In fact, without separating 
the contributions of the positive parity baryon states from the negative parity baryon 
states explicitly, the two criteria (pole dominance and convergence of the operator prod- 
uct expansion) for choosing the Borel parameter T 2 and threshold parameter sq in the 
conventional QCD sum rules do not work efficiently; we maybe (or maybe not) choose the 
Borel windows where the contaminations from the negative parity doubly heavy baryon 
states are large. With suitable Borel parameters, we can choose the tensor structures ji or 
1 to study the masses and pole residues freely. If we choose the tensor structure 70 + 1, 
the contaminations from the negative parity doubly heavy baryon states are excluded 
explicitly. 

Those doubly heavy baryon states E cc , f2 cc , and fl^b maybe observed at the Teva- 
tron, the LHCb and the PANDA, especially at the LHCb. For example, the 5 CC and O cc can 
be produced at the high energy pp or pp collisions through the gluon-gluon fusion mecha- 
nism and the intrinsic charm mechanisms, g+g — > (cc)[ 3 Sijg+c+c, g+g — > (cc)[ So]6+c+c, 
g + c -)• (cc)[ 3 5i]3 + c, g + c -)• (cc)^ S ]q + c 1 c + c ^ (cc)[ 3 5i]§ + g, c + c -)• (cc)^ S ] e + g, 
where the (cc)[ 3 5i]3 (in color anti-triplet 3) and (cc)[ 1 5o]6 (in color sextet 6) are two 
possible 5-wave configurations of the doubly charmed diquark pair (cc) inside the baryon 
states H cc and O cc [371 ES EHl 113 E] . The LHCb is a dedicated b and c-physics precision 
experiment at the LHC (large hadron collider). The LHC will be the world's most copi- 
ous source of the b hadrons, and a complete spectrum of the b hadrons will be available 
through gluon fusion. Furthermore, once reasonable values of the pole residues An and A3 
are obtained, we can take them as basic input parameters and study the revelent hadronic 
processes with the QCD sum rules. 

4 Conclusion 

In this article, we study the ^ + doubly heavy baryon states Oqq and Sqq by subtracting 
the contributions from the corresponding i doubly heavy baryon states with the QCD 
sum rules, and make reasonable predictions for their masses. Those doubly heavy baryon 
states maybe observed at the Tevatron, the LHCb and the PANDA, especially at the 
LHCb. Once reasonable values of the pole residues A^ and A~ are obtained, we can take 
them as basic input parameters and study the revelent hadronic processes with the QCD 
sum rules. 
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